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1. Introduction 

One of the most important problems of the modern string theory is a search for new 
nontrivial backgrounds and the study of the underlying string dynamics. In most of the 
cases the target space metric of such backgrounds is curved and often it incorporates 
the black hole singularities. In the super string theories, the supersymmetry allows for 
some interesting nontrivial solutions which are stable and exact. But the string theory on 
such backgrounds is usually an extremely complicated sigma-model, very difficult even to 
formulate it explicitly, not to mention studying quantitatively its dynamics. 

The two-dimensional bosonic string theory is a rare case of sigma-model where such 
a dynamics is integrable, at least for some particular backgrounds, including the dilatonic 
black hole background. A physically transparent way to study the perturbative (one loop) 
string theory around such a background is provided by the CFT approach. However, once 
we want to achieve more ambitious quantitative goals, especially in analyzing higher loops 
or multipoint correlators, we have to address ourselves to the matrix model approach to the 
2D string theory proposed in [3j] in the form of the matrix quantum mechanics (MQM) of 
. The string theory has been constructed as the collective field theory (the Das- Jevicki- 
Sakita theory), in which the only excitation, the massless "tachyon", was related to the 
eigenvalue density of the matrix field. 

In the framework of MQM it is difficult to formulate directly a string theory in a 
nontrivial background metric since the operators perturbing the metric do not have a 
simple realization. However, we can perturb the theory by other operators, like tachyon 
or winding operators. We hope that such a perturbation can also produce a curved metric 
but in an indirect way. An example of such a duality was given in || where the 2D black 
hole background is induced by a winding mode perturbation. 

The 2D string theory in the simplest, translational invariant background (the linear 
dilaton background) is described by the singlet sector of MQM with inversed quadratic 
potential. In the singlet sector of the Hilbert space, the eigenvalues have Fermi statistics, 
and the problem of calculating the S'-matrix of the string theory tachyon becomes a rather 
standard problem in a quantum theory of a one-dimensional nonrelativistic free fermionic 
field. The tree- level S'-matrix can be extracted by considering the propagation of "pulses" 
along the Fermi surface and their reflection off the "Liouville wall" 0,[5|. 

The very existence of a formulation in terms of free fermions indicates that the 2D 
string theory should be integrable. For example, the fermionic picture allows to calculate 
the exact S'-matrix elements. Each S'-matrix element can be associated with a single 
fermionic loop with a number of external lines [Q . One can then expect that the theory is 
also solvable in a nontrivial, time-dependent background generated by a finite tachyonic 
source. Dijkgraaf, Moore and Plesser |7j demonstrated that this is indeed the case when 
the allowed momenta form a lattice as in the case of the compactified Euclidean theory. 
In it has been shown that the string theory compactified at any radius R possesses the 
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integrable structure of the Toda lattice hierarchy. However, this important observation 
had not been, until recently, really exploited. The Toda structure is too general, and it 
becomes really of use only if accompanied by an initial condition or a constraint (a " string 
equation" ), which eliminates all the solutions but one. Thus the first results for a nontrivial 
background in case of a general radius were obtained as a perturbative expansion in the 
tachyon source ||]. 

Recently, the Toda integrable structure of the compactified Euclidean 2D string theory 
was rediscovered [p|,pT0H and used to find the explicit solution of the theory [jl^JO] • Finally, 



the string equation at an arbitrary compactification radius was found in ||12|| . These pa- 
pers studied the T-dual formulation of the string theory, where instead of the discrete 
spectrum of tachyon excitations the winding modes of the string around the compactifi- 
cation circle (Berezinski-Kosterlitz-Thouless vortices) were used. These modes appear in 
the non-singlet sectors of MQM [113 1 and can be generated by imposing twisted periodic 
boundary conditions |T4[ . 



In this paper we return to the study of the Toda integrable structure of tachyon 
excitations of 2D string theory originated in 0. It describes special perturbations of the 
ground state within the singlet sector of the MQM. We will construct the Lax operators 
as operators in the Hilbert space of the singlet sector of MQM. We will be able to find 
an interpretation of the Toda spectral parameters in terms of the coordinates of the two- 
dimensional target phase space and interpret the solutions of the Toda hierarchy in terms 
of the shape of the Fermi sea. In particular, we find the explicit shape of the Fermi sea for 
the Sine-Liouville string theory. 

We can give two different interpretations of our problem. The first one is that of a 
2D string theory in Minkowski space in a non-stationary background. The simplest, time- 
independent ground state of the theory is characterized by a condensation of the constant 
tachyon mode, which is controlled by the cosmological constant /i. Here we will study 
more general, time-dependent backgrounds characterized by a set of coupling constants 
t± n associated with non-trivial tachyon modes with purely imaginary energies E n = in/R, 
where R is a real number. Since the incoming/outgoing tachyons with imaginary energies 
have wave functions exponentially decreasing/increasing with time, such a ground state 
will contain only incoming tachyons in the far past and only outgoing tachyons in the far 
future. In other words, the right and left vacua are replaced by coherent states depending 
correspondingly on the couplings t n and t- n (n = 1,2, ...), which will be identified with 
the "times" of the Toda hierarchy. The coherent states of bosons modify the asymptotics 
of the profile of the Fermi sea at far past and future, without changing the number of 
fermions. The flow of the fermionic liquid is no more stationary, but its time dependence 
happens to be quite trivial, and the profiles of the Fermi sea at different moments of time 
are related by Lorentz boosts. 
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The second interpretation of our analysis is that we consider perturbations of an 
Euclidean 2D string theory in which the Euclidean time is compactified at some length 
(3 = 2nR. Together with the cosmological term, we allow all possible vertex operators 
with momenta n/ R, n G 2. The simplest case of such a perturbation is the Sine-Liouville 
theory. This case will be considered in detail and the shape of the Fermi sea produced by 
the Sine-Liouville perturbation will be found explicitly. The solution we have found exhibits 
interesting thermo dynamical properties, which may be relevant to the thermodynamics of 
the string theory on the dilatonic 2D black hole background. 

Let us mention also another possible application of our analysis: the fermionic system 
that appears in MQM is similar to the problem of two-dimensional fermions in a strong 
transverse magnetic field [^5j. The electrons filling the first Landau level form stationary 
droplets of Fermi liquid, similar to that of the eigenvalues in the phase space of MQM. The 
form of such droplets is also described by the Toda hierarchy [^6[. Our problem might be 
related to the situation in which two such droplets are about to touch, which can happen 
at a saddle point of the effective potential []T7| . 

The paper is organized as follows. In section 2 we will remind the CFT formulation 
of the 2D string theory. In section 3 the MQM in the inverted harmonic potential in the 
"light-cone" phase space variables is formulated and the free fermion structure of its singlet 
sector is revealed. In section 4 the one particle wave functions are studied. In section 5, 
after the description of the fermionic ground state, tachyonic perturbations are introduced 
and the equations defining the corresponding time-dependent profile of the Fermi sea are 
derived. In section 6 the Toda integrable structure of the perturbations restricted to a 
lattice of Euclidean momenta is derived directly from the Schroedinger picture for the free 
fermions. In section 7 we recover the solution of the Sine-Liouville model and describe 
explicitly its Fermi sea profile. In section 8 we reproduce the free energy of the perturbed 
system from the profile of the Fermi sea. The section 9 is devoted to conclusions and in 
section 10 we discuss some problems and future perspectives. In particular, we propose a 
3-matrix model describing the 2D string theory perturbed by both tachyon and winding 
modes. In the appendix the one particle wave functions of the type II model, defined on 
both sides of the quadratic potential, in the "light cone" formalism are presented. 

2. Tachyon and winding modes in 2D Euclidean string theory 

The 2D string theory is defined by Polyakov Euclidean string action 



where the bosonic field x(a) describes the embedding of the string into the Euclidean time 
dimension and g ab is a world sheet metric. In the conformal gauge g ab = e^ a ^g ab , where 




(2.1) 
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g ab is a background metric, the dilaton field cf> becomes dynamical due to the conformal 
anomaly and the world-sheet CFT action takes the familiar Liouville form 



S = ^-[ d 2 a [(9x) 2 + (d(j)) 2 -m<t> + ne~ 2 ^ + ghosts]. (2.2) 

J world sheet 

This action describes the unperturbed linear dilaton string background corresponding to 
the flat 2D target space (x, <$>). 

In the target space this theory possesses only one propagating degree of freedom which 
corresponds to tachyon field. If the Euclidean time is compactified to a circle of radius R, 
x(a) = x(a) + 2ttR, the spectrum of admissible momenta is discrete: p n = n/R, n E Z. 
In this case there is also a discrete spectrum of vortex operators, describing the winding 
modes (Kosterlitz-Thouless vortices) on the world sheet. A vortex of charge q m = mR 
corresponds to a discontinuity 2-KmR of the time variable around a point on the world 
sheet .0 The explicit expressions of the vertex operators V p and the vortex operators V q in 
terms of the position field x = xr + xl and its dual x = xr — xl are 

v _ r (N) f -i P x (\p\-2)4> 

P ~n-\P\)J ' (23) 

Any background of the compactified 2D string theory can be obtained (at least in the 
case of irrational R) by a perturbation of the action (|2.2j ) with both vertex and vortex 
operators 

S = S Q + Y,{t n V n + t n V n ). (2.4) 

Such a theory possesses the T-duality symmetry of the non-perturbed theory x <-> x, 
R — > 1/-R, fx — > fx/R ]T^] if one also exchanges the couplings as t n <-> t n . Another general 
feature of the theory ( |2.4f) is the existence of the physical scaling of various couplings, 
including the string coupling g s , with respect to the cosmo logical coupling fi. It can be 
found from the zero mode shift of the dilaton <p — > 4> + 0o- I n this way we obtain 

t„~A* 1 " 4H/Jl . h"!* 1 -* 11 ™, 9s -t*- 1 . (2.5) 

The scaling ( |2.5| ) allows to conclude that at the compactification radius lower than the 
Berezinski-Kosterlitz-Thouless radius Rrt = 1/2 all vertex operators are irrelevant. In 



1 At rational values of R there exist additional physical operators (similar to discrete states 
at the self-dual radius R = 1) containing the derivatives of fields x{a) and 4>(u), which we will 
ignore in this paper. 
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the interval 1/2 < R < 1 the only relevant momenta are p = ±1/R. The theory perturbed 
by such operators looks as the Sine-Gordon theory coupled to 2d gravity ( "Sine-Liouville" 
theory): 



- / d 2 ( [(8x) 2 + (8<P) 2 - 41l<f) + ixe- 2(t> + Xe^' 2 ^ cos(x/i2)]. (2.6) 
in J 



It was conjectured by [[| that at R = 2/3 and \i = the "Sine-Liouville" is dual to the 



SL(2,C) 



coset model with central charge c = ~ 1 = 26, which describes the 2D 



_3k 

SU{2)xU(l) 

string theory in the black hole ("cigar") background. 

If we go to the Minkowski space, the perturbation (|2.4| ) is made by tachyons with 
purely imaginary momenta. In this case there are more general perturbations, generated 
by tachyons with any real energies. In the next section we will introduce the matrix 
formulation of the string theory in Minkowski space using light-cone coordinates. 



3. Matrix Quantum Mechanics in light cone formulation 

The 2D string theory in Minkowski time appears as the collective field theory for 
the large- N limit of MQM in the inverted oscillator potential Jl Jl9 , 20| , 21 j22 ] . The matrix 
Hamiltonian is 

H= -Tr(P 2 -M 2 ), (3.1) 
2 

where P = —id/dM and My is an iV x N hermitian matrix variable. The cosmological 
constant \i in (|2.2|) is introduced as a "chemical potential" coupled to the size of the 
matrix N, which should be considered as a dynamical variable. In this formulation the 
time coordinate of the string target space coincides with the MQM time (or its Euclidean 
analogue x = —it in the compactified version of the MQM) and that of the Liouville field 
is related to the spectral variable of the random matrix. 

The tachyon modes of the string theory are represented by the asymptotic states of 
the collective theory. The scattering operators with real energy E and describing left- and 
right-moving waves, respectively, are given by 

Tf = e~ lEt Tr (M + P) iE , T E = e~ iEt Tr (M — P)~ iE . (3.2) 

These operators can be used to construct the in- and out-states of scattering theory. 
Namely, for an in-state, one needs a left-moving wave while the out-state is necessarily 
given by a right-moving one. The vertex operators ( |2.3| ) are tachyons with purely imagi- 
nary energies and are therefore represented by the following chiral operators in the MQM 



J e - pt Tr(M + P)H, p>o, 
p ~* \e" pt Tr(M-P)IPl, p<0. 



(3.3) 
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Since we are interested in perturbations with the chiral operators (|3.3| ), it is natural to 
perform a canonical transformation to light-cone variables 



M±P , - 

Z ± = — (3.4) 

and write the matrix Hamiltonian as 

H = -±Tr(Z + Z-+Z-Z+), (3.5) 

where the matrix operators Z± obey the canonical commutation relation 

[(£+)}, (Z_)f\ = -i 5\5*. (3.6) 

Define the right and left Hilbert spaces as the spaces of functions of Z+ and Z_ only, with 
the scalar product 

($±|$' ± ) = J d N2 Z ± <f>±(Z±)& ± (Z ± ). (3.7) 

The operator of coordinate in the right Hilbert space is the momentum operator in the 
left one and the wave functions in the Z+ and Z_ representations are related by a Fourier 
transform. The second-order Schrodinger equation associated with the Hamiltonian ( |3.1] ) 
becomes a first-order one when written in the i-representations 

&t*±(Z ± ,t) = ^Tr \ Z ± -^- + ^pj^ ± (Z ± ,t), (3.8) 

whose general solution is 

$±(Z±,t) = e^ N2t ^ { l\e^Z ± ). (3.9) 

The Hilbert space decomposes into a direct sum of eigenspaces labeled by the irre- 
ducible representations of SU(N), which are invariant under the action of the Hamiltonian 
( |3.5|) . The functions 

$(r) = ) belonging to given irreducible representation r 

transform as 

^(ntz±fi) = J2 D uW$£ J \z ± ), (3.io) 
j 

(r) 

where Djj is the group matrix element in representation r and /, J are the representation 
indices. The wave functions transforming according to given irreducible representation 
depend only on the iV eigenvalues z ± 1 , . . . , z ± N and the Hamiltonian ( |3.5| ) reduces to its 
radial part 

v ^ N 
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A potential advantage of the light-cone approach is that the Hamiltonian ( |3.11|) does not 
contain any angular part, which is not the case for the standard one ( |3.1| ), whose angular 
term induces a Calogero-like interaction. 

In the scalar product (|3.7|) , the angular part can also be integrated out, leaving only 



the trace over the representation indices: 

= J2 I n^ ± , fc A 2 (, ± )^ J) (,J^ J) (, ± ), (3.12) 

J k 

where A(z ± ) is the Vandermonde determinant. If we define 

^\z ± ) = A(z ± )^\z ± ), (3.13) 

these determinants disappear from the scalar product. The Hamiltonian for the new 
functions &±\z ± ) takes the same form as ( |3.11| ) but with a different constant term: 

#o = X>±. At,* +1 /2). (3.14) 

k 

In the singlet representation, the wave function ty± (z ± ) = \]/2 mglet ' ) (z ± ) is a completely 
antisymmetric scalar function. The scalar product in the singlet representation is given by 

k 

Thus the singlet sector describes a system of iV free fermions. The singlet eigenfunctions 
of the Hamiltonian ( |3.14| ) are represented by Slater determinants of one-particle eigenfunc- 



tions discussed in the next section. 

It is known ]T3|,|T4|,|TU[ that unlike singlet, which is free of vortices, the adjoint repre- 



sentation contains string states with a vortex-antivortex pair, and higher representations 
contain higher number of such pairs. In what follows we will concentrate on the fermionic 
system describing the singlet sector of the matrix model. We will start from the proper- 
ties of the ground state of the model, representing the unperturbed 2D string background 
and then go over to the perturbed fermionic states describing the (time-dependent) back- 
grounds characterized by nonzero expectation values of some vertex operators. 



4. Eigenfunctions and fermionic scattering 

To study the system of non-interacting fermions we have to start with one-particle 
eigenfunctions. The one-particle Hamiltonian in the light-cone variables of the previous 
section can be written as 

2 

where z, turn out to be canonically conjugated variables 



H = ~-{z + z_ +z_z + ), (4.1) 



[z +1 z_] = -i. (4.2) 

We can work either in z + or in z_ representation, where the theory is defined in terms of 
fermionic fields ip ± (z ± ) respectively. General solutions of the Schrodinger equation with 
the Hamiltonian ( f4.1| ) written in these representations take the form 

^ ± (z ± ,t)=z- 1 / 2 f ± (e^ z± ). (4.3) 

There are two versions of the theory, referred in || as theories of type I and II. In 
the theory of type I the eigenvalues A& of the original matrix field M are restricted to 
the positive half-line. The theory of type II is defined on the whole real axis A. Such a 
string theory splits, at the perturbative level, into two disconnected string theories of type 
I. Here we will consider the fermion eigenfunctions in the theory of type I. The theory of 
type II will be considered in Appendix A. 

In the light cone formalism it is natural to replace the restriction A > with z ± > 0, 
which again does not affect the perturbative behavior. In this case the solutions with a 
given energy are tp^(z ± ,t) = e~ lEt i/j^(z ± ) with 

^(s ± ) = 4=s? B -*- (4-4) 

The functions ( f4.4j ) with E real form two complete systems of 5-function normalized or- 
thonormal states 



POO 

{n\^l')= / dz ± J^)^{z ± )=5{E-E') 1 (4.5) 
Jo 



dE ^(z ± )^(z' ± ) = 8(z ± -z' + ). (4.6) 



As any two representations of a quantum mechanical system, the z + and z_ represen- 
tations are related by a unitary operator S, which in our case is nothing but the Fourier 
transformation on the half-line. The latter is defined by the integral 



[£ty + KO= / dz + K(z_,z + )^ + (z + ), (4.7) 



o 
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where there are two choices for the kernel: 



K(z_,z + ) = )J-cos{z_z + ) or K(z_, z + ) = iy — sin(z_z + ). (4-8) 

The sine and the cosine kernels describe two possible theories, which differ only on the 
non-perturbative level.! Let us choose the cosine kernel in (^j8|) and evaluate the action of 
the S'-operator on the eigenfunctions. This integral is essentially the defining integral for 
the T-function: 

1 r°° 

[^VfKag = -J dz ± cos{z + z_)zf E ~* = K(±E)^(z^, (4.9) 

where 

K(E) = v/^cosh (^(V 2 - E )) + V 2 )- (4-10) 

(The sine kernel would give the same result, but with cosh replaced by sinh). The factor 
7Z(E) is a pure phase 

k(e)h(e) = n{-E)n{E) = 1, (4.11) 

which proves the unitarity of the operator S. The operator S relates the incoming and 
the outgoing waves and therefore can be interpreted as the fermionic scattering matrix. 
The factor 1Z(E) is identical to the the reflection coefficient (the "bounce factor" of |§), 
characterizing the scattering off the upside-down oscillator potential. In the standard 
scattering picture, the reflection coefficient is extracted by comparing the incoming and 
outgoing waves in the large-A asymptotics of the exact eigenfunction of the inversed oscil- 
lator Hamiltonian H = —^{d\ + A 2 ). 

It follows from the above discussion that the scattering amplitude between two arbi- 
trary in and out states is given by the integral with the Fourier kernel 



r<*> (4.12) 

(ip_\K\ip + ) = / dz + dz_ip_(z_) K(z_,z + )ip + (z + ). 



The integral (|4.12|) can be interpreted as a scalar product between the in and out states. 
By ( |4.5| ) and ( fO| ) one finds that the in and out eigenfunctions satisfy the orthogonality 
relation 

(ip^\K\iJj®'} = 1Z(E)8(E — E 1 ). (4.13) 



2 The fact that there are two choices for the kernel can be explained as follows. In order 
to define the theory of type I for the original second-order Hamiltonian (3.1), we should fix the 
boundary condition at A = 0, and there are two linearly independent boundary conditions. 
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5. String theory backgrounds as profiles of the Fermi sea 

The ground state of the MQM is obtained by filling all energy levels up to some 
fixed Fermi energy which we choose to be Ep = —fx. Quasiclassically every energy level 
corresponds to a certain trajectory in the phase space of z + , z_ variables. The Fermi sea 
can be viewed as a stack of all classical trajectories with E < Ep and the ground state 
is completely characterized by the curve representing the trajectory of the fermion with 
highest energy Ep. For the Hamiltonian ( |4. 1| ) all trajectories are hyperboles z + z_ = —E 
and the profile of the Fermi sea is given by 

z + z_ = fx. (5-1) 

In the quasiclassical limit the phase-space density of fermions is either or 1. Then 
the low lying collective excitations are represented by deformations of the Fermi surface, 
i.e., is the boundary of the region in the phase space filled by fermions. At any moment 
of time such deformation can be obtained by replacing the constant (i on the right hand 
side of (|5.1|) with a function of z + and z_ 

z + z_ = M(z + ,z_). (5.2) 

In contrast to the ground state which is stationary, an excited state given by a generic 
function M leads to a time dependent profile. However, this dependence is quite trivial: 
since the Fermi surface consists of free fermions each moving according to its classical 
trajectory z ± (t) = e ±t z ± (0), where z ± (0) is the initial data, we can always replace ( |5T^ ) 
by the equation for the initial shape. So the evolution of a shape in time is simply its 
homogeneous extension with the factor e* along the z + axis and a homogeneous contraction 
with the same factor along the z_ axis. 

Below we will find equations that determine the shape of the Fermi sea for a generic 
perturbation with tachyon operators. Our analysis is in the spirit of the Polchinski's deriva- 
tion of the tree- level S'-matrix [|J, but we will consider finite, and not only infinitesimal 
perturbation. 

In terms of the Fermi liquid, the incoming and the outgoing states are defined by the 
asymptotics of the profile of the Fermi sea at z + ^> z_ and z_ ^> z + , correspondingly. If 
we want to consider such a perturbation as a new fermion ground state, we should change 
the one-fermion wave functions. The new wave functions are related to the old ones by a 
phase factor 

Vl(z ± ) = e T ^ ±(z ± ;S V|(^ ± ), (5-3) 

whose asymptotics at large z ± characterizes the incoming/outgoing tachyon state. We 
split the phase into three terms 

<p±(z ± ;E) = V ± (z ± ) + ^(E)+v ± (z ±] E), (5.4) 
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where the potentials V± are fixed smooth functions vanishing at z ± = 0, while the term 
v± vanishing at infinity and the constant <p are to be determined. Thus, the potentials V± 
fix unambiguously the perturbation. The time evolution of these states with the original 
Hamiltonian ( |4.1|) is determined by the eq. ([Of). To see that the state flOp contains 
incomming and outgoing tachyons, it is enough to note that it arises as a coherent state 
of vertex operators ( |3.3|) acting on the unperturbed wave function (|4.4| ). 



Since the functions and V?^ should describe the same one-fermion state, the Fourier 
transform ([Op should be diagonal in their basis. We fix the zero mode (p so that 

Stf£ = ^. (5.5) 

This condition can also be expressed as the orthonormality of in and out eigenfunctions 

(^ E _-\K\^l + ) = 6{E+-E-), (5.6) 

with respect to the scalar product ( f4.12| ). This requirement fixes the exact form of the 
wave functions. Let us look at this problem in the quasiclassical limit E± — > oo. The 



scalar product in (|5.6|) is written as 



oo 



(*M< + > = ^= / d ^^ C os{z + z_)e-^ {z + ) -^ z - ) z^ E -zf+. (5.7) 



IZ , z 







At the quasiclassical level it can be evaluated by the saddle point approximation. One 
obtains two equations for the saddle pointi 

z + z_ = —E± + z ± d(p±(z ± ). (5.8) 



Generically, the two equations (§^) define two different curves in the z + -z_ plane, and 
their compatibility can render at most a discrete set of saddle points. However, if the two 
solutions define the same curve, we obtain a whole line of saddle points, which implies the 
existence of a zero mode in the double integral ( |5.7| ). The contribution of the zero mode 
is infinite, amounting to the ^-functional orthogonality relation. Thus, the orthonormality 
condition is reduced to the compatibility of two equations ( p.8| ) at E + = E_. 

For example, in the absence of perturbations (ip± = 0), the saddle point equations 
z + z_ = —E± are inconsistent, unless E+ = E-. For equal energies they coincide leading 
to a zero mode in the double integral ( [4. 13| ) . The resulting saddle point equation gives the 
classical hyperbolic trajectory in the phase space of an individual fermion. If E+ ^ E-, 
the integrand of (|5.7|) is a rapidly oscillating function and the integral is zero. 



3 We assume that only one exponent of cos(z + z_) gives a contribution. The result will be 
shortly justified from another point of view. 
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The equations ( |5.8|) and the requirement of their compatibility can be obtained from 
another point of view. The perturbed wave function (|5.3|) can be interpreted as an eigen- 
state of a new, perturbed Hamiltonian. Indeed, the functions ( |5.3| ) are not eigenstates of 
the original Hamiltonian Q4.1| ) in ±-representations Hq = ^i(z ± d z± + 1/2), but for each 
given E they are evidently eigenstates of the operators 

HZ=H£ + z ± d<p ± {z ± ;E), (5.9) 

where <p± contain so far unknown functions v±(z ± ; E). However, the operators H± de- 
pend on the energy E through these functions. Therefore, they can not be considered as 
Hamiltonians. But one can define the Hamiltonians as solutions of the equations 

H± = H± + z ± d<p ± (z ± ;H ± ). (5.10) 

Then all functions ( |5.3| ) are their eigenstates. 

The orthonormality condition ( |5.6|) can be equivalently rewritten as the condition 
that the Hamiltonians H± define the action of the same self-adjoint operator H in the 
i-representations. In the quasiclassical limit, this is equivalent to the coincidence of the 
phase space trajectories associated with H + and f/_: 

H+(z + ,z_) = E H-(z + ,z_) = E. (5.11) 

This condition is equivalent to the compatibility of two equations (|5.8|) . 



Note that with respect to the time r defined by the new Hamiltonian, the time de- 
pendence of the states characterized by the wave functions ( |5.3|) is given by e~ iEr . This 
corresponds to a stationary flow of the Fermi liquid. The profile of the Fermi sea coincides 
with the classical trajectory of the fermion with the highest energy Ep = Its equation 
can be written in two forms similar to ( |5.2| ) which should be consistent 

z + z_ = M±(z ± ) = fi + z ± d<p±(z ± ;-fj). (5.12) 

These equations are the non-compact version of the equations arising in the conformal map 
problem which is a semiclassical description of compact Fermi droplet |16[]. In that case 



the potential must be an entire function which is not required in our case. The functional 
equation ( 5.12 ) contains all the information of the tachyon interactions in the tree level 



string theory. To proceed further with our analysis, we should concretize the form of 
the perturbing potentials V±. We will show in the next section that the perturbations 
produced by vertex operators are integrable. 
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6. Integrable flows associated with vertex operators 

6.1. Lax formalism, Toda Lattice structure and string equations 

Now we restrict ourselves to the time dependent coherent states made of tachyons 
with discrete Euclidean momenta p n = nj R. This spectrum of momenta arises when the 
system compactified on a circle of length 2tcR, or heated to the temperature T = 1/(2ttR). 
Such perturbations are described by the potentials of the following form 



In this section we will show that such a deformation is exactly solvable, being generated by a 
system of commuting flows H n associated with the coupling constants t± n . The associated 
integrable structure is that of a constrained Toda Lattice hierarchy. The method is very 
similar to the standard Lax formalism of Toda theory, but we will not assume that the 
reader is familiar with this subject. It is based on rewriting all operators in the energy 
representation, in which the Fourier transformation S is diagonal as we required in the 
previous section. The energy E will play the role of the coordinate along the lattice 
formed by the allowed energies E n = ip n of the tachyons. 

Let us start with the operators z ± = z ± and <9 T = ^j-, which are related as 



V±(z ± ) = RY,t± k z 



± 



k/R 



(6.1) 



fc>i 



— id. 



Sz.S 



id. 



'+ = 



S- X z S. 



(6.2) 



The Heisenberg commutation relation [d±, z ± ] = 1 leads to the operator identity 



z_Sz + S 1 — Sz + S 1 z_=i. (6-3) 

Further, the Hamiltonian (fOJ) in the z ± -representation = ^fi(z ± d z± +1/2) is expressed 
in terms of z, and S as 




It follows from the identities 



z ± ^(z ± ) 



E^(z ± ), 



(6.5) 



that the operators Hq, S and z, are represented in the basis of the non-perturbed wave 



functions (4.4) by 



H^ = E, z ± =u ±1 , S = K(E)l = Kl, (6.6) 
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where Lip^_ = ip^ and cD is the shift operator 

u = e~ idE . (6.7) 

In order to have a closed operator algebra, we should also define the action of the operators 
S ±x on the functions i(}± Tni (z ± ) = z^il)±(z ± ). To satisfy the identities (|6.2[), we should 
define S ,±1 by the cosine kernel for n even and by the sine kernel for n odd. Therefore 
Sip^^) = VJ^\z ± ) where K' = 1Z'(E) is given by ( pHl) with cosh replaced by sinh. 
The algebraic relations ( |6.3| ) and ( |6.4j) are transformed into algebraic relations among 



the operators E, £j and 1Z in the .E-spacel 

K^ujU'Cj- 1 = -E + i/2, 

, . , - . (6.8) 
u TZ ujTZ = —E — i/2. 

This is equivalent, according to the remark above, to the functional constraint 

K{E) = -{i/2 + E)K'{E + i), (6.9) 

which is evidently satisfied. To simplify the further discussion, in the rest of this sec- 
tion we will identify the functions 1Z(E) and TZ'(E), thus neglecting the non-perturbative 
terms ~ e nE . Note however that all statements made below can be proved without this 
identification. 

Our aim is to extend the i?-space representation for the basis of wave functions (|4.4| ) 
perturbed by a phase factor e ±lLp± as in eqs. (|5.3|) -(|5^4T), with V± given by flO]). We 
assume that the phase ip± can be expanded, for sufficiently large Laurent series 

fc>i fc>i 

This assumption will be justified by the subsequent analysis. With the special form fl6.10j ) 
of the perturbing phase, there exist unitary "dressing operators" W± acting in the .E-space, 
which transform the "bare" wave functions into the "dressed" ones 

= e T<v±(*±)^| = W±^|. (6.11) 

It is evident from ( p.lUp and ( |6.6D that the dressing operators can be expressed as series 
in the shift operator ui with independent coefficients 

W± = e^/ 2 (l + X>±*a T * /J M e^^*^**'* (6.12) 

\ k>l I 



Note that we should inverse the order of the operators. 
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In the basis of the perturbed wave functions (|6.11| ), the operators of canonical coor- 
dinates z ± are represented by two Lax operators 

L± = W± u ±l W± 1 = e T^/2 c2 ,±i e ±#/2 j\ + £ a±fc ^=Ffc/-R j (6.13) 

and the Hamiltonians are represented (up to a change of sign) by the so-called Orlov- 
Shulman operators 



M± = -Wi^Wi 1 = ^ kt ±k L k J R -E + J2 v± k L ± h/R . (6.14) 

fc>i fc>i 



In the above formulas all coefficients are functions of E and t = {t± k }kLi- The Lax and 
Orlov-Shulman operators act on the wave functions as 

L±*U*±) = *±*±{z±), (6-15) 

M±*£(z ± ) = I kt± k z k J R + v±kZ~ k/R I *J(z ± ) (6.16) 

\k>l k>l J 

and satisfy the Lax equations 

[L±,M±} = ±iL ± , (6.17) 
which are the dressed version of the relation [cu ±:L , — E] = ±iu . 



The structure of constrained Toda hierarchy follows from the requirement (|5.5| ) , which 
means that the action of the S operator on the perturbed wave functions is totally com- 
pensated by the dressing operators: 

W- = W+TZ. (6.18) 

The condition ( |6.18| ) defines both the Toda structure and a constraint which plays the role 
of an initial condition for the PDE of the Toda hierarchy. 

The Toda structure implies that the tachyon operators generating the perturbation 
are represented in the -E-space by an infinite set of commuting flows. To show this, we 
evaluate the variations of the Lax operators with respect to the coupling constants t n . 
From the definition ( |6.13| ) we have 



dt n L± = [H n ,L±], (6.19) 
where the operators H n are related to the dressing operators as 

H n = (d tn W+)W+ l = (d tn W-)WZ l . (6.20) 
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The two representations of the flows H n are equivalent by virtue of the relation ( |6.18|) . 
It is important that 7Z does not depend on t n 's. A more explicit expression in terms 
of the Lax operators is derived by the following standard argument. Let us consider 
the case n > 0. From the explicit form of the dressing operators ( |6.12|) it is clear that 
H n = W+ u) n / R W+ 1 + negative powers of cD 1,/i? . The variation of t n will change only the 
coefficients of the expansions (|6.13|) of the Lax operators, preserving their general form. But 



it is clear that if the expansion of H n contained negative powers of u) 1 /^, its commutator 
with L- would create extra powers ui~ 1 ~ k / R . Therefore 

H ±n = (L n ± /R )> + ±(L n ± /R ) , n>0, (6.21) 

where the symbol ( )> means the positive (negative) parts of the series in the shift 
operator Cj and ( )o means the constant part. By a similar argument one shows that the 
Lax equations ( |6.19| ) are equivalent to the zero-curvature conditions 

d tm H n - d tn H m - [H m , H n ] = 0. (6.22) 



The equations (|6.19|) and (|6.21|) ensure that the perturbations related with the couplings 



t n are described by the Toda integrable structure. 

The Toda structure leads to an infinite set of PDE's for the coefficients w n of the 
dressing operators, the first of which is the Toda equation for the zero mode of the dressing 
operators 

2 JL -^—(j) = e^W-^iE+i/R) - e i^{E-i/R)-i${E) _ 23 ) 

The uniqueness of the solution is assured by appropriate boundary conditions or additional 
constraints, known also as string equations. In our case the string equation can be obtained 
by dressing the equation ( |6.8| ) using the formula ( |6.18| ), which leads to the following relation 
between the Lax and Orlov-Shulman operators 

L+L-=M+ + i/2, L-L+=M--i/2. (6.24) 

Similarly, the identity 1ZE = E1Z implies 

M_ = M+. (6.25) 

Thus, the perturbations of the one-fermion wave functions of the form (|6.11|) are 



described by a constrained Toda lattice hierarchy. In the T-dual formulation the same 



result was proved recently using the standard definition of the Toda Lax operators ||12|| . 
The standard Lax and Orlov-Shulman operators are related to L± and M± by 

L = L+ /i? , L = L l J R , M = M + , M = M_. (6.26) 
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This operators satisfy .R-dependent constraints 
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L R L R = M + i/2, L H L R = M -i/2, M = M. 



(6.27) 



(The string equations for integer values of R have been discussed previously in [23.24,25,2^].) 

The integrability allows to find the unknown coefficients in the Laurent expansion 
of the phase ( |6.10| ). All the information about them is contained in the so called r- 
function, which is the generating function for the vector fields corresponding to the flows 
H n . Its existence follows from the zero-curvature condition (|6.22| ) and can be proved in 



the standard way ||27|1 . The coefficients are related to the r-function as follows 



<91ogr 
dt n 



4>{E) = i\og 



r(E + i/2R) 
r(E-i/2R) 



(6.28) 



The r-function can thus be considered as the generating function for the perturbations by 
vertex operators. As we will show in sec. 8, the logarithm of the r-function gives the free 
energy of the Euclidean theory compactified at radius R. 

6.2. The dispersionless (quasiclassical) limit 

Let us consider the quasiclassical limit E — ► — oo when all allowed energies are large 
and negative. In this limit the integrable structure described above reduces to the disper- 
sionless Toda hierarchy |]28,^9|,27l, where the operators E and Gj can be considered as a 
pair of classical canonical variables with Poisson bracket 



{cu, E} = u. 



(6.29) 



Similarly, all operators become c-functions of these variables. The Lax operators can be 
identified, by eq. ( |6.15| ), with the phase space coordinates z ± . The functions z ± (u,E) are 
given by eq. 



(6.30) 



fc>i 

where we have introduced the "string susceptibility" 

X = -Rd E 4> = d 2 E \ogr. 
The equations ( |6.24j) and ( |6.25|) in the quasiclassical limit give 

{z_,z + } = 1, 
z + z_ = M + = M_. 



(6.31) 

(6.32) 
(6.33) 
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It follows from the first equation, (|6.32| ), that eq. (|6.30 ) defines a canonical transforma- 
tion relating the phase-space coordinates z + ,z_ to logo; and E. Moreover, the dressing 
procedure itself can be interpreted as a canonical transformation between z + , z_ and the 
"bare" coordinates 

uj + = ^f-E w, uj- = V^E/u. (6.34) 
The second equation, ( |6.33j ), which is a deformation of the relation 

u+u- = -E (6.35) 

satisfied by the "bare" variables, yields the form of the perturbed one-fermion trajectories. 

Let us consider the highest trajectory that defines the shape of the Fermi sea. The 
variable E then should be interpreted as the Fermi energy, which is equal to (minus) the 
chemical potential, Ep = —fi. Taking into account the expression for M± through z ± given 
by the r.h.s. of ( |6.16| ), we can expand the right-hand side either in z + or in z_ 

z + z^=J2 kt±k z k J R + + $>±fc Z ± k/R - ( 6 - 36 ) 
fc>i fc>i 

(These expansions are of course convergent only for sufficiently large values of z ± .) Eq. 
( |6.36|) , which is the dispersionless string equation, is a particular case of the relation ( |5.12|) , 
which was derived for an arbitrary potential. In this way we found a direct geometrical 
interpretation for the dispersionless string equation in terms of the profile of the Fermi sea. 

A simple procedure to calculate the coefficients a±k has been suggested in |l~2| . First, 
note that if all t±k with k > n vanish, the sum in ( |6.30| ) can be restricted to k < n. Then 
it is enough to substitute the expressions ( J6.3U[ ), with E = —fi, in the profile equations 



( |6.36| ) and compare the coefficients in front of uj ±k ^ R . The result will give the generating 
function of one point correlators, which is contained in the inverse function u(z ± ). In the 
next section we will apply this method to calculate the profile of the Fermi sea in the 
important case of the Sine-Liouville string theory. 



7. Solution for Sine-Gordon coupled to gravity 

The simplest nontrivial string theory with time-dependent background is the Sine- 
Gordon theory coupled to gravity known also as Sine-Liouville theory. It is obtained by 
perturbing with the lowest couplings t\ and t—\. It can be easily seen that in this case the 
expansion ( |6.30| ) consists of only two terms 

z± = e"^ x w ±1 (l + a±w T *). (7.1) 
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The procedure described in []T2| gives the following result for the susceptibility x an d the 
coefficients a±: 



1 \ 2B-1 ., 

l__)t 1 t_ ie ^x = l5 



a± = tzpi e h 2 x . 



(7.2) 



The first equation in ( |7.2| ) for the susceptibility X = log r was found in ggl and it was 
shown to reproduce the Moore's expansion of the free energy F = logr The equation 
( [7.1| ) with a± given in ( |7.2| ) was first found in [ITJ in the form 



e^± - x ±e h± = 1, 



(7.3) 



where /i±(z ± ) is the generating function of the one-point correlators 



It is related to iv(z ± ) as follows 



d 2 F 



n=l 



u>{z ± ) = e 2RA -z ± e 



-Rh±(z ± ) 



(7.4) 



(7.5) 



The equation (|7.2|) for the susceptibility imposes some restrictions on the allowed 
values of the parameters t±\ and /i. Let us consider the interval of the radii \ < R < 1, 
which is the most interesting from the point of view of the conjectured duality between the 
Sine-Liouville theory and the 2D string theory in the black hole background at R = 2/3 
0. In this interval, if we choose tit-i < 0, a real solution for x exists only for fx > 0. 
On the other hand, if > 0, the allowed interval is \x > (j, c , where the critical value is 
negative: 



1 

2 

R 



1 

R- 1 



{tit. 



(7.6) 



Let us note that ( |7.6|) was interpreted in |31J as a critical point of the pure 2D gravity 
type. 

The solution ( |7.1] )-( [fT2|) allows us to find the explicit shape of the Fermi sea in the 
phase space. In the model of the type I considered so far there is an additional restriction 
for the admissible values of parameters since the boundary of the Fermi sea should be 
entirely in the positive quadrant z ± > 0. This forces us to take t±i > 0. The general 
situation is shown on Fig. la. The unperturbed profile corresponds to the hyperbola 
asymptotically approaching the z ± axis, whereas the perturbed curves deviate from the 



Since the papers [|H||rj|l!,|3rj considered the case of vortex perturbations, the duality trans- 
formation R — > 1/R should be performed before comparing the formulas. 
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Fig. 1: Profiles of the Fermi ± p) in the theory of the type I at 

R = 2/3. Fig. la contains several profiles corresponding to t\ = £_i = 2 and 
values of \x starting from \x c = — 1 with a step 40. For comparison, we also 
drew the unperturbed profile for \i = 100. Fig. lb shows three moments of 
the time evolution of the critical profile at \i = — 1. 

axes by a power law. We see that there is a critical value of \i 1 where the contour forms 
a spike. It coincides with the critical point \i = [i c given by (|7.6| ). At this point the 
quasiclassical description breaks down. On the Fig. lb, the physical time evolution of a 
profile is demonstrated. 

In the same way one can find the solution in the classical limit of the theory of the 
type II described in Appendix A. In this case we can introduce two pairs of perturbing 
potentials describing the asymptotics of the wave functions at z ± — > oo and z ± — > — oo. 
For sufficiently large \i the Fermi sea consists of two connected components and the theory 
decomposes into two theories of type I. However, in contrast to the previous case, there are 
no restrictions on the signs of the coupling constants. When \i decreases, the two Fermi 
seas merge together at some critical value \i* . This leads to interesting (for example, from 
the point of view of the Hall effect) phenomena, which we intend to discuss elsewhere. 
Here we will only mention that, depending on the choice of couplings, it can happen that 
for some interval of /i around the point \i* , the Toda description is not applicable. 

8. Free energy of the perturbed background 

As we mentioned, the perturbations considered above appear in a theory at the finite 
temperature T = 1/(3, (3 = 2ttR. The free energy per unit volume T is related to the 
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partition function Z by 

Z = e~ pT . (8.1) 

The free energy is considered as a function of the chemical potential //, so that the number 
of fermions is given by 

N = dF/dft. (8.2) 

We have seen that any perturbation can be characterized by the profile of the Fermi 
sea. For a generic profile, the number of fermions is given by the volume of the domain in 
the phase space occupied by the Fermi liquid 

N = -|- J J dz + dz_. (8.3) 

Fermi sea 

It is implied that the integral is regularized by introducing a cut-off at distance a/A. For 
the unperturbed ground state ( |5.1| ), one reproduces from (|8.2| ) and ( |8.3|) the well known 
result for the (universal part of the) free energy 

T=^ 2 \og^ (8.4) 



The authors of |7| identify the partition function ( |8.1| ) for the tachyon backgrounds 
studied in sect. 6 with the r-function of the Toda hierarchy, defined by ( |6.28| ). In this 
section we will reproduce this statement by direct integration over the Fermi sea, adjusting 
the above derivation of ( |8.4j ) to the case of a general profile of the form ( |6.36| ). 

Since the change of variables described by eq. ( |6.30|) is canonical (see eq. (|6.32| )), one 
can rewrite the integral (|8.3| ) as 

The limits of integration over u are determined by the cut-off.li If we put it at the distance 
z ± = y/A, they can be found from the equations 



z ± (u; ± (E),E) = VA. (8.6) 
Taking the derivative with respect to fi, we obtain 

1 uj , (—it) 



6 The integral over E is also bounded from below, but we do not need to specify the boundary 
explicitly. 
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It is enough to keep only the leading order in the cut-off A for the boundary values oj . 
From ( |8.6| ) and ( |6.30| ) we find to this order 

zb 1 /2 

u ± = (ke^ R ) . (8.8) 



Combining (|8.7|) and (|8.8|) , we find 

1 . , 1 



d ll N = -—\ogk-— X . (8.9) 
Ztv ZtvK 



Taking into account the relation (8.2) between the free energy and N, we see that 



jr = _I logT5 (3 = 2ttR. (8.10) 

P 

As a result, if we compactify the theory at the time circle of the length f3, our free energy 
coincides with the logarithm of the To da r-function. 

The perturbed flow of the Fermi liquid is non-stationary with respect to the physical 
time t associated with the original Hamiltonian Hq. In sec. 5 we constructed a new 
Hamiltonian, which preserves the shape of the Fermi sea. The latter coincides with the 
phase-space trajectory ( |5.11|) at E = —fx. Using the fact that d fl H(z + , z_) = 0, one can 



check that 

F=(H(z + ,z_)+n). (8.11) 

Note that the function H(z + , z_) is obtained by solving the profile equations ( |6.36| ) with 
respect to \i. One can see, looking at the large z ± asymptotics of the profile equations, 
that it can be written in the following compact form: 

H{z + ,z_) =H + J2 k (tkH k + t-kH- k ). (8.12) 
fc>i 

Here the Hamiltonians H n (z + , z_ ; t) are given by equations ( |6.21| ), with u and E expressed 
as functions of z + and z_ through eq. (|6.30| ). 



9. Conclusions 

In this paper we studied the 2D string theory in the presence of arbitrarily strong 
tachyonic perturbations. In the matrix model language, we dealt with the singlet sector of 
the MQM described by a system of free fermions in inverted quadratic potential. In the 
quasiclassical limit the state of the system is described by the shape of the domain in the 
phase space occupied by the fermionic liquid (the profile of the Fermi sea). In this limit 
we formulated functional equations, which determine the shape of the Fermi sea, given its 
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asymptotics at two infinities. There is some analogy between our problem and the problem 



of conformal maps studied in [^ , |3^ , |34| , |35| . 

This allows to calculate various physical quantities related to the perturbed shape, 
such as the free energy of the compactified Euclidean theory, in terms of the parameters 
defining the asymptotic shapes. 

For a particular case of perturbations generated by vertex operators (in- and out-going 
with discrete equally spaced imaginary energies), the system solves the Toda lattice hier- 
archy where the Toda "times" correspond to the couplings characterizing the asymptotics 
of the profile of the Fermi sea. Using the "light cone" representation in the phase space, 
we developed the Lax formalism of the constrained Toda hierarchy, where the constraint 
coincides with the string equation found recently in [[12|| . In the dispersionless limit, the 
Lax formalism reproduces the functional equations for the shape of the perturbed Fermi 
sea and allows to give a direct physical meaning to the Toda mathematical quantities. 
In particular, we identified the two spectral parameters of the Toda system with the two 
chiral coordinates z ± in the target phase space. 

Let us list the other possible applications of our formalism: 

• Calculation of the tachyon scattering matrix for 2D string theory in nontrivial back- 
grounds. 

• Analysis of the systems with more complicated than Sine-Liouville backgrounds. In 
particular, the quasiclassical analysis of the sections 4 and 7 are not limited to any lattice of 
tachyon charges and can be used for a mixture of any commensurate or non-commensurate 
charges. It provides a way of studying the whole space of possible tachyonic backgrounds 
of the 2D string theory. 

• Investigation of the thermodynamics of the 2D string theory in Sine-Liouville type tachy- 
onic backgrounds characterized by a lattice of discrete Matsubara energies |30] , |36] . This is 
supposed to be the thermodynamics of the dilatonic black hole, according to the conjecture 
of §• 

• The type 2 theory can be used to model the situation where either two quantum Hall 
droplets approach each other or one droplet is about to split in two. 



10. Problems and proposals 

We list also some unsolved problems, which could be approached by our formalism: 

• To find the metric of the target space for the Sine-Liouville string theory. 

• To construct the discrete states of the 2D string theory (at any rational compactification 
radius) in the framework of our formalism. 

• A few important unsolved problems concern the relation between the MQM and CFT 
formulations of the 2D string theory. We still don't have a precise mapping of the vertex 
and vortex operators of the MQM to the corresponding operators of CFT, in spite of the 
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useful suggestions of the papers |37], Q. As the result, we did not manage to match 
the Sine-Liouville/Black Hole correlators calculated by |§, [j38l , and []39j from the CFT 
approach with the corresponding correlators found in |11|] from the Toda approach to the 
MQM. Another related question is how can the target space integrable structure be seen 
in the CFT formulation of 2D string theory. 

• Description and quantitative analysis of the mixed (vertex&vortex) perturbations in the 
2D string theory by means of the MQM. We hope that the "light-cone" formalism might 
be appropriate to study such perturbations, in spite of the absence of integrability.0 

Finally, keeping in mind this interesting problem, let us propose here a 3-matrix model 
whose grand canonical partition function Z(fi,t n ,t n ) = Y^N = o e ~ 2nRfJ,N ^N(tn^n) is the 
generating function of correlators for both types of perturbations: the parameters t n are 
the couplings of the vortex operators, in the same way as £ n 's are the couplings of the 
tachyon (vertex) operators throughout this paper. As usual, the chemical potential // 
plays the role of the string coupling. The partition function at fixed iV is represented by 
the integral over two hermitian matrices Z± and one unitary matrix, the holonomy factor 
O around the circle 



Z N (t n ,t n )= [VZ + ] {tn} [VZ_] {t _ n} [dn] {i±n} e 



Tr (Z-Z + -qZ_Q,Z + Q,- r 



(10.1) 



where we denoted q = e l2 ' KR and introduced the following matrix integration measures 



[dZ ± ] {t±n} = dZ ± e*£»>o 



[dn] a 



{t± n } 



[dn] 



SU(N) 



— J 1 



t„ Tr Q r 



(10.2) 



10.3) 



This 3-matrix integral is nothing but the Euclidean partition function of the upside-down 
matrix oscillator on a circle of length 2nR with twisted periodic boundary conditions 
Z±(P) = & Z±(0)Q, needed for the introduction of the winding modes (vortices) in the 
MQM model of the 2D string theory [p~4 j . It is trivial to show that in particular case of 
all t n = this 3-matrix model reduces to the model with only vortex perturbations of the 
It can be also shown 13F 



paper 
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that in the case of all i n = it reduces to the Euclidean 
version of the model with only tachyonic perturbations studied in this paper. The crucial 
point here is that the sources for tachyon perturbations are introduced in a single point on 
the time circle, similarly to the perturbation imposing the initial conditions on the wave 
function ( [Blf ). The model ( |10.1|) reduces the problem of the study of the backgrounds 



7 An interesting proposal to describe nontrivial string backgrounds on the CFT side was given 
long ago in [4C] and elaborated in [41],[42|. It suggested to parameterize the moduli space of the 
2D string theory perturbed in both sectors, by a conifold. This structure may be hidden in the 3 
matrix model proposed below. 
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of the compactified 2D string theory with arbitrary tachyon and winding sources, to a 
three-matrix integral. 

We learned from R. Dijkgraaf and C. Vafa that they discovered a similar Toda struc- 
ture in the c = 1 type theory arising in connection to the topological strings. 
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Appendix A. Theory of type II 

In the theory of type II, the fermions are defined on the whole real line. In this case 
one should introduce two sets of functions describing, in the quasiclassical limit, fermions 
at different sides of the potential. They are defined on right (left) semi-axis by 



±iE- 



i 



1 2 

V>f (Q = -!= f ± (z+>o), 



i (-O 



±iE- 



(A.l) 



±,<V ±! V2^ Vl + e 2 * E ± 



and the continuation to the other semi-axis is performed according to the rule z ± — > e T7rt ; 
This gives, for z ± > 



^ ±> {-z ± )=±ie E ^l > {z ± ) 1 



(A.2) 



The functions ( |A.1|) satisfy the orthonormality and completeness conditions 

(iPl a \^ b ) = 6 ab 8(E-E f ), (A.3) 

/DO 

dE YMM)KS Z ' ± ) = S ( z ± - <=). (A.4) 

where the indices a, b = (>, <) label the right and left wave functions and the scalar 

product is defined as in ( |4.5| ) but with the integral over the whole real axis. 

The unitary operator relating z ± representations is given by the Fourier kernel on the 

whole line K(z_,z.) = —k=e lz + z - 
< ' v ^ 

poo 

[S^ + ](z_) = / dz + K(z_,z + )i; + (z + ). (A.5) 
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It acts on the eigenfunctions ( |A.1|) by the following matrix 

^VfJ^HE^WVf^), S(E)=K(E)(_. 1 enE "'f*), (A.6) 
where 

K{E) = -^e-^ E ~^ 2) T{iE + 1/2). (A.7) 

V 27T 

In this case the reflection coefficient 7Z(E) is not a pure phase because of the tunneling 
through the potential described by the off-diagonal elements of the S'-matrix. Nevertheless, 
the whole S'-matrix is unitary 

S^(E)S(E) = 1. (A.8) 
The scalar product between left and right states is given by 

/oo 
dz + dz_ JJT)K{z_ ,z + )iP + {z + ). ( A.9) 

-oo 

It is clear that the matrix of scalar products of states ( |A.1| ) coincides with S{E)5{E — E'). 
The corresponding completeness condition is 

dE 5>- a (0 {S- l ) ab r + , b {z + ) = -r=e~ iz + z -. (A.10) 

a,6 V2n 
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